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ABSTRACT 

We study the non-linear evolution of a dust ellipsoid, embedded in a Fried- 
mann fiat background universe, in order to determine the evolution of the density 
of the ellipsoid as the perturbation to it related detaches from general expansion 
and begins to collapse. We show that while the growth rate of the density contrast 
of a mass element is enhanced by the shear in agreement with Hoffman (1986a), 
the angular momentum acquired by the ellipsoid has the right magnitude to 
counterbalance the effect of the shear. The result confirms the previrialization 
conjecture (Peebles & Groth 1976; Davis & Peebles 1977; Peebles 1990) by show- 
ing that initial asphericities and tidal interactions begin to slow down the collapse 
after the system has broken away from the general expansion. 

Subject headings: cosmology: theory — galaxies: formation 

1. Introduction 

While cosmologists generally believe that structures in the universe grew by gravita- 
tional instability from smaller inhomogeneities present at the epoch of decoupling, there is 
disagreement on several details of the model. One of these is the role of asphericity in the 
collapse of perturbations and structure formation. 

According to the previrialization conjecture (Peebles & Groth 1976, Davis & Peebles 1977, 
Peebles 1990), initial asphericities and tidal interactions between neighboring density fluc- 
tuations induce significant non-radial motions which oppose the collapse. This means that 
virialized clumps form later, with respect to the predictions of the linear perturbation theory 
or the spherical collapse model, and that the initial density contrast, needed to obtain a given 
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final density contrast, must be larger than that for an isolated spherical fiuctuation. This 
kind of conclusion was supported by Barrow & Silk (1981), Szalay & Silk (1983), Villumsen 
& Davis (1986), Bond & Myers (1993a,b) and Lokas et al. (1996). 

In particular Barrow & Silk (1981) and Szalay & Silk (1983) pointed out that non-radial 
motions would slow the rate of growth of the density contrast by lowering the peculiar ve- 
locity and suppress collapse once the system detaches from general expansion. Villumsen 
& Davis (1986) gave examples of the growth of non-radial motions in N-body simulations. 
Arguments based on a numerical least-action method lead Peebles (1990) to the conclusion 
that irregularities in the mass distribution, together with external tides, induce non-radial 
motions that slow down the collapse. Lokas et al. (1996) used N-body simulations and 
a weakly non-linear perturbative approach to study previrialization. They concluded that 
when the slope of the initial power spectrum is n > —1, non-linear tidal interactions slow 
down the growth of density fluctuations and the magnitude of the effect increases when n is 
increased. 

Opposite conclusions were obtained by Hoffman (1986a, 1989), Evrard & Crone (1992), 
Bertschinger & Jain (1994) and Monaco (1995). In particular Hoffman (1986a,1989), using 
the quasi-hnear (QL) approximation (Zel'dovich 1970; Zel'dovich & Novikov 1983) showed 
that the shear affects the dynamics of collapsing objects and it leads to infall velocities that 
are larger than in the case of non-shearing ones. Bertschinger & Jain (1994) put this result 
in theorem form, according to which spherical perturbations are the slowest in collapsing. 
The N-body simulations by Evrard & Crone (1992) did not reproduce previrialization effect, 
but the reason is due to the fact that they assumed an n = — 1 spectrum, differently from the 
n = one used by Peebles (1990) that reproduced the effect. If n < —1 the peculiar grav- 
itational acceleration, g oc it!~("+^)/^, diverges at large R and the gravitational acceleration 
moves the fluid more or less uniformly, generating bulk flows rather than shearing motions. 
Therefore, its collapse will be similar to that of an isolated spherical clump. If n > —1, the 
dominant sources of acceleration are local, small-scale inhomogeneities and tidal effects will 
tend to generate non-radial motions and resist gravitational collapse. 

In a more recent paper. Audit et al. (1997) have proposed some analytic prescriptions 
to compute the collapse time along the second and the third principal axes of an ellipsoid, 

by means of the 'fuzzy' threshold approach. They point out that the formation of real 
virialized clumps must correspond to the third axis collapse and that the collapse along this 
axis is slowed down by the effect of the shear rather than be accelerated by it, in contrast 
to its effect on the first axis collapse. They conclude that spherical collapse is the fastest, in 
disagreement with Bertschinger & Jain's theorem. This result is in agreement with Peebles 
(1990). 

In this paper, we address this controversy by following the evolution of a dust ellipsoid 
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in an expanding universe. We shall use a model of Nariai & Pujimoto (1972), that makes 
possible to study separately the effect of the shear, E, and that of angular momentum, L, 
on the protostructure evolution. 

The paper is organized as follows. In section 2, we describe the model used; in section 
3 we calculate the angular momentum of the ellipsoid at an intermediate time between the 
turn-around of the first and third axis, and in section 4 we describe the parameters and 
initial conditions used. Section 5 is devoted to the discussion of results and section 6 to 
conclusions. 



2. Ellipsoid model for the collapse 

In order to determine the evolution of the density in an ellipsoid of ideal fluid at zero 
pressure, wc follow Nariai & Fujimoto (1972) and Barrow &; Silk (1981). In what follows, 
we shall study the evolution of the density, p, of an ellipsoid embedded in a pressureless 
background cosmology with zero curvature characterized by a background density, pb, and 
expansion parameter a{t): 

As showed by Nariai & Fujimoto (1972), performing two transformations of coordinates from 
the co-moving frame {x'^} to the inertial frame x^ — {t, x') with: 

x' = a(t)x' (2) 

and then from the last to a non-inertial system of reference, {x''^}, rotating with angular 
velocity fli with respect to the inertial frame {x^}, the Newtonian hydrodynamical equations 
of continuity, motion and Poisson are: 
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where p and p are the density and pressure perturbations and (henceforth designating x[ and 
Vl with Xi and Vi, respectively): 
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Considering a rotating ellipsoid having uniform density, p — p{t), the velocity field is 
given by: 

Vi = aijXj (7) 

with: 

1 /«! 6:2 0:3 



(8) 

where the shear tensor, Sy, and the vorticity vector, a;k, are respectively given by: 



77 \ \ 5ij + Sij + £ijka;k 

3 \ai a2 ^3 



_ 1 1 

^ij = 2 ("ij + ~ 3«kk5ij (9) 

a;i = SijkQijk (10) 

The term aj represent the i — th principal semi-axis of the ellipsoid. Assuming that the 
rotation velocity possesses only an 2:3 component and that the initial vorticity (note that 
here the term 'initial vorticity' shall not be interpreted as primordial vorticity, which is zero 
before orbit crossing, but as the vorticity acquired after shell-crossing) has no components 
in the directions of xi and X2, then the equations of motion for the principal axes of the 
ellipsoid and that for the evolution of density are: 



ai = — pb 1 - -aia2a3Ui ai - -GMUiai + -3 (11) 

6 \ Z J I ^q;^ + 0,2) 



OL2 = — Pb 1 - -Ot\OL20i-iU2 Oi2 - -GMU2a2 + -3 (12) 

6 \ Z ) I (^Q,^ + Q,2j 

Al'kG / 3 \ 3 

cig = — I 1 - -axOL20L3U3 ) tta - -GMU^a^ (13) 



g-^(g)'-4.Gp-S^+ ,, ^0 (14) 

p o \p J aia2{a.i + a2Y 

(Nariai & Fujimoto 1972; Barrow & Silk 1981), where L is the angular momentum of the 
ellipsoid 



V^(x) = nK + ^) (16) 

i=l 
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M — —paia2a2 — constant (17) 
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and 



3 



= = E? + E^ + E^ 



(18) 
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(see equation (3.22) of Nariai & Fujimoto (1972)). Equation (14), even if not strictly neces- 
sary to describe the evolution of the ellipsoid, is very useful because it allows us to qualita- 
tively understand how the gravitational instability process is modified by the rotation and 
shear anisotropy. If E = L = 0, we are reconducted to the spherically symmetric case of no 
rotation. In this case the density reaches a maximum value of 97r^/16pb and after the system 
recoUapses. The shear, E^, acts in the same sense of gravity making collapse easier, while the 
angular momentum L acts in the opposite sense to self-gravity, Gp, making it easier to resist 
gravitational collapse. Obviously equation (14) could be also used to obtain the evolution 
of the density after calculating E^, and L^. As shown by Nariai & Fujimoto (1972), E^ is 
given by equation (18) and E; can be obtained by means of equation (9) and equation (3.4) 
of Nariai & Fujimoto (1972) as follows: 



This last result shows even clearly that equation (11, 12, 13) form a closed system of equations 
giving the evolution of the ellipsoid and the shear. 

Then the evolution of the density can be obtained (once the angular momentum is 
known) in two ways: 

1) by integrating equation (11, 12, 13) to get the evolution of the semi-axes. Then E^ can 
be calculated through equation (19) and equation (18) and finally the density evolution is 
obtained by integrating equation (14). 

2) By integrating equation (11, 12, 13) to get the evolution of the semi-axes and then using 
equation (17). 

It is useful to remark that while the procedure 2) is simpler than 1) in the case E 7^ 0, in 
the case L = E = 0, it is simpler to use the procedure 1) because, in this case, we a priori 
know that E = and consequently we have only to perform the integration of equation (14). 
If, otherwise, we wanted to use the procedure indicated in the point 2) we should integrate 
equation (11, 12, 13) and also impose the condition that E^ = using equation (19), with 
a consequent complication of the calculations. In any case, in the paper, we performed the 
calculations following both the procedures, in order to check the consistency of the results. 



O CCi CCj CKk 



(19) 



A fundamental point to remark is the following: the Nariai & Fujimoto (1972) equations 
give a description of the evolution and collapse of an ellipsoid only if the ellipsoid has acquired 
somehow angular momentum. For the reasons described in the following, we use this model 
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to study the evolution of the eUipsoid from the epoch of turn-around of the first axis on. 
An eUipsoid can have angular momentum for two different reasons: 

1) the axes of the ellipsoid and that of the shear of the velocity field have an appropriate 
misalignement, or in other terms the principal axes of the inertia tensor arc not aligned with 
the principal axes of the deformation tensor (see Catelan & Theuns 1996a,b). In this case, 
the ellipsoid can have angular momentum even if vorticity is zero. We are not interested in 
this particular case. 

2) the system has a non-zero vorticity. Unfortunately, we now that according to Kelvin's 
theorem, if the initial velocity field is irrotational, i.e. curl-free, then, it should remain 
irrotational also in the nonlinear regime. However, since the collapse of a protostructure 
is a violent phenomenon, the conditions of Kelvin's circulation theorem should be violated 
(Chernin 1970). Then, there are two possibility for vorticity generation (see Sasaki & Kasai 
1997): 

a) acquisition of vorticity by the formation of shock fronts in the protostructure (pan- 
cake), in correspondence of shell- crossing (Doroshkevich 1970). Analytical studies by Pichon 
& Bernardeau (1999) have also shown that vorticity generation becomes significant at the 
scales 3 — 4h~^Mpc, and increases with decreasing scale; 

b) acquisition of angular momentum by means of the tidal torques (Hoylc, 1949; Peebles 
1969; Hoffman 1986b; Ryden & Gunn 1987). Current analytical description of vorticity and 
spin growth by tidal torques turn out to depend on a free parameter, i.e. the time when tidal 
torques cut off. This parameter has been related to cither the beginning of the decoupling 
from the Hubble flow (5 ~ 1) or the turn-around epoch (time when expansion halts) (for 
the spherical collapse model) (Andriani & Caimmi 1994). Numerical simulations (Barnes 
& Efstathiou 1987) have shown that after decouphng from the Hubble flow there are no 
substantial increment in angular momentum. 

Summarizing, we know (and assume) that from the linear phase to shell-crossing the 
vorticity is zero, the ellipsoid has no rotation. The perturbation is subject to the gravitational 
fleld of matter inside the ellipsoid, which tends to make it collapse to a pancake, and to the 
tidal fleld of the matter outside, which cancels the effects of the interior gravitational fleld. 
As a result, the ellipsoid expands with the rest of the universe and preserves its shape until 
it enters the nonlinear phase (Barrow & Silk 1981; White & Silk 1979). When it reaches a 
density contrast 5 ^ 1 it detaches from Hubble expansion, and the distribution of matter 
of the ellipsoid tends to develope nonradial motions (Peebles 1980), (even if the axes of the 
ellipsoid and that of the shear of the velocity fleld have not an appropriate misalignement). 

To take account of the rotation acquired, we identify the flnal angular momentum of 
the ellipsoid with that acquired at the maximum of expansion of the object (Peebles 1969; 
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Catelan & Theuns 1996a,b). 

This assumption is justified by the fact that after the maximum expansion time the 
angular momentum stops growing, becoming less sensitive to tidal couplings (Peebles 1969; 
Barnes & Efstathiou 1987). 

In other words, we assume that the total angular momentum is acquired before the 
ellipsoid collapse, since the tidal torques are much less effective afterwards (Peebles 1969; 
Catelan & Theuns 1996a,b). While for a spherical perturbation this time is well defined (it 
is the turn-around epoch) , for an aspherical perturbation this epoch is not well defined and 
the system should be followed until the long axis turns around (Eisenstein & Loeb 1995), 
since the acquisition of angular momentum is important until the collapse of this axis. To 
simplify things, we choose a mean value of time, tyi, between the turnaround of the shortest 
axis and that of the longest (see Hoffman 1986b) 

3. Anguleir momentum at tu 

As previously remarked, if we want that equation (11, 12, 13) constitute a closed system 
of equations to get the ellipsoid evolution, we need the angular momentum L. As stressed 
previously, we need only the value of the angular momentum of the virialized structure, that 
can be well approximated (see the above discussion) by its value at the time tu- 

The effect of tidal torques on structures evolution has been studied in several papers 
especially in connection with the origin of galaxies rotation. The explanation of galaxies 
spins gain through tidal torques was pioneered by Hoyle (1949) in the context of a collapsing 
protogalaxy. Peebles (1969) considered the process in the context of an expanding world 
model, showing that the angular momentum gained by the matter in a random comoving 
Eulerian sphere grows at the second order in proportion to t^/^ (in a Einstein-de Sitter uni- 
verse), since the proto-galaxy was still a small perturbation, while in the non- linear stage 
the growth rate of an oblate homogeneous spheroid decreases with time as t~^. 
White (1984) considered an analysis by Doroshkevich (1970) that showed that the angular 
momentum of galaxies grows to first order in proportion to t and that Peebles's result is 
a consequence of the spherical symmetry imposed to the model. White showed that the 
angular momentum of a Lagrangian sphere does not grow either in the first or in the second 
order, while the angular momentum of a non-spherical volume grows to the first order in 
agreement with Doroshkevich's result. 

Another way to study the acquisition of angular momentum by a proto-structure is that 
followed by Ryden (1988) (hereafter R88) and Eisenstein & Loeb (1995). Following Eisen- 
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stein & Loeb (1995), we separate the universe into two disjoint parts: the collapsing region, 
characterized by having high density, and the rest of the universe. The boundary between 
these two regions is taken to be a sphere centered on the origin. As usual, in the follow- 
ing, we denote with p(x), being x the position vector, the density as function of space and 
S{x.) = ^Mz£b_ fjj^g gravitational force exerted on the spherical central region by the ex- 
ternal universe can be calculated by expanding the potential, $(x), in spherical harmonics. 
Assuming that the sphere has radius R, we have: 



oo , 

4:71 



*(^) = E ^j—i E ^Ux)Y,^{e, ct>)x' (20) 
where Fim arc spherical harmonics and the tidal moments, aim, are given by: 

poo 

aira{x)^fh> / Y,^{eA)p{^)s-'-^dh (21) 

JR 

In this approach the proto-structure is divided into a series of mass shells and the torque 
on each mass shell is computed separately. The density profile of each proto-structure is 
approximated by the superposition of a spherical profile, 5(r), and the same Gaussian density 
field which is found far from the peak, £(r), which provides the quadrupole moment of the 
proto-structure. To the first order, the asphericity about a given peak can be represented 
writing the initial density in the form: 

p(r) = pb[l + 5(r)][l+£(r)] (22) 

(Ryden & Gunn 1987; R88; Peebles 1980, equation (18.5)) where e{v) satisfies the following 
equations: 

(k.r) = P{k) (23) 

being P{k) the power spectrum, and: 

(£(r)) = 0, {e{vf) = m (24) 

(Ryden & Gunn 1987; Peebles 1980), where <> indicates a mean value of the physical quan- 
tity considered and ^ the two-point correlation function (note that the previous equations 
are obtained in the lowest order approximation. The limits of the same are described in 
Ryden & Gunn 1987). The torque on a thin spherical shell of internal radius x is given by: 

r{x) = ^ / £(x)xxv$(x)dQ (25) 

where Mgh — 47rpb [1 + ^{x)] x'^Sx. Before going on, I want to recall that we are interested 
in the acquisition of angular momentum from the inner region, and for this purpose we take 
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account only of the I — 2 (quadrupole) term. In fact, the I — term produces no force, while 
the dipole (/ = 1) cannot change the shape or induce any rotation of the inner region. As 
shown by Eisenstein & Loeb (1995), in the standard CDM scenario the dipole is generated 
at large scales, so the object we are studying and its neighborhood move as bulk flow with 
the consequence that the angular distribution of matter will be very small, then the dipole 
terms can be ignored. Because of the isotropy of the random field, £(x). Equation (25) can 
be written as: 

< kr >= a/(30)^ [< a2m{xf >< q2m{xf >-< a2m{x)q2m{x) >^Y^^ (26) 

As stressed in the next section, following Eisenstein & Loeb (1995), the integration of the 
equations of motion shall be ended at some time before the inner external tidal shell (i.e., 
the innermost shell of the part of the universe outside the sphere containing the ellipsoid) 
collapses. Then the inner region behaves as a density peak. This last point is an important 
one in the development of the present paper. 

An important question to ask, before going on, regards the role of triaxiality of the el- 
lipsoid (density peak) in generating a quadrupole moment. Equation (26) takes into account 
the quadrupole moment coming from the secondary perturbation near the peak. The density 
distribution around the inner region is characterized (see Equation (22)) by a mean spher- 
ical distribution, S, and a random isotropic field. In reality the central region is a triaxial 
ellipsoid. It is then important to evaluate the contribution to the quadrupole moment due 
to the triaxiality. Remembering that the quadrupole moments are given by: 

= / Y*^{e, <p)s'p{s)dh = [ Y*^{e, <p)e{^)dn (27) 

J\r\<R 47r J 

and approximating the density profile as: 

5(x) =< d{x) > Spherical +iyf{x)A{e,p) (28) 

being < 5{x) > spherical the mean spherical profile, — ^ the peak height and a the r.m.s. 
value of S. The function A{e,p) of the triaxiality parameters, e and p, is given by: 

A(e, p) = 3e(l - sin^ 9 - sin^ 9 sin^ 0) + p(l - 3 sin^ 9 cos^ 0) (29) 

while the function f{x) is given (R88) by: 

where ^, a and i?* are respectively the two-point correlation function, the mass variance and 
a parameter connected to the spectral moments (see Bardeen et al. 1986, equation (4.6d), 
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hereafter BBKS). Substituting equation (28) and equation (29) in equation (27) it is easy to 
show that the sum of the mean quadrupole moments due to triaxiahty is: 

^ ^ < q^^ix) >= vx^f{x) (^i-v^(e -p) + l-^/^{Ze + p)^ (31) 

which must be compared with that produced by the secondary perturbations, e: 

< q2^{xf >= ^^Ml I k'P{k)Ukxfdk (32) 

where j2 is the Bessel function of order 2. The values of e and p can be obtained from the 
distribution of elhpticity and prolateness (BBKS, equation (7.6) and figure 7) or for i/ > 2 
by: 

e = T7^ (33) 

v^x [1 + 6/(5x2)]^/' 

and 

p 2 (34) 

5x4 [1 + 6/(5x2)]2 V J 

(BBKS equation (7.7)), where x is given in BBKS (equation (6.13)). In the case of a peak 
with u = 3, wc have e ^ 0.15, p ~ 0.014 while for peaks having i> = 2 and v = 1 they are 
respectively given by e ~ 0.2, p ~ 0.03 and e ~ 0.25 p ~ 0.04. 

As shown in figure 1, for a 3(T profile, the source of quadrupole moment due to triaxiahty 
is less important than that produced by the random perturbations e in all the proto-structure, 
except in the central regions where the quadrupole moment due to triaxiahty is comparable 
in magnitude to that due to secondary perturbations. In other words, the triaxiahty has 
a significant effect only in the very central regions, which contains no more than a few 
percent of the total mass and where the acquisition of angular momentum is negligible. It 
follows that the triaxiahty can be ignored while computing both expansion and spin growth 
(R88). Moreover, as observed by Eisenstein & Loeb (1995), the ellipsoid model does better 
in describing low shear regions (having higher values of z/), whose collapse is more spherical 
and then the effects of triaxiahty are less evident. Just this peaks, having at least z/ > 2, 
shall be studied in this paper. In any case, even if the triaxiahty was not negligible it should 
contribute to increment the acquisition of angular momentum (Eisenstein & Loeb 1995), and 
finally to a larger effect on the density evolution, (i.e., a larger reduction of the growing rate 
of the density) . 

In order to find the total angular momentum imparted to a mass shell by tidal torques, it 
is necessary to know the time dependence of the torque. This can be done connecting g2m and 
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a2m to parameters of the spherical collapse model (Eisenstein & Loeb 1995 (equation (32), 
R88 (equation (32) and (34)). Following R88 we have: 

and 

a2mid) = a2m,o 77 So{9 - sm9) 3 (36) 



3/ 



The collapse parameter 9 is given by: 



t{9) = -to6o''^\e-sm9) (37) 
Equation (35) and (36), by means of equation (26), give to us the tidal torque: 

3 3 (^-«i^^)^'^'Vi(^)-(|)/2(^) 

where fi{9) and f2{d) arc given in R88 (Eq. 31), Tq and 6q = arc respectively the torque 
and the mean fractional density excess inside the shell, as measured at current epoch tg- The 
angular momentum acquired during expansion can then be obtained integrating the torque 
over time: 

£ = /rW> (39) 

As remarked in the previous section, the angular momentum obtained from equation(39) 
is evaluated at the time t^- Then the calculation of the angular momentum can be solved 
by means of equation (39), once we have made a choose for the power spectrum. With the 
power spectrum and the parameters given in the next section and ior a u — 2 peak, the 
model gives a value of 2.5 x lO^^gcm^/s. 

Since the calculation of the angular momentum is fundamental for the evolution of the 
ellipsoid, it is worthwhile to comment on the validity of the calculation and the result. 

To start with, we want to recall that, independently from the he calculation followed 
in order to get the angular momentum, we need only its final value. Then it is important 
to compare the result obtained with our approach with those obtained following different 
approaches. An interesting comparison is that with the result of Catelan & Theuns (1996a). 
In their paper, they calculated the angular momentum at maximum expansion time (see 
their equations (31)-(32)) and compared it with previous theoretical and observational es- 
timates. Assuming the same value of mass u used to obtain our previously quoted result 
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(2.5 X lO^^gcmVs) and the same value of redshift {z — 3) and distribution of the final an- 
gular momentum /f adopted by Catelan & Theuns (1996a), we get a value for the angular 
momentum of 2 x lO^^gcm^/s. This last result is in good agreement with ours and is well 
in line with previous theoretical estimates (Peebles 1969; Heavens & Peacock 1988) and 
numerical simulations (Fall 1983). It is obvious that the approach of this paper or that of 
Catelan & Theuns (1996a) cannot predict the very final stages of evolution when clumps 
merge and interact non-linearly and in addition dissipative processes may play an important 
role as well. In any case, the value of the final angular momentum, as obtained from the 
extrapolation of the hnear theory, is typically a factor of ~ 3 larger than the final spin of the 
non-hnear object (Barnes & Efstathiou 1987; Prenk 1987). The effects of this discrepancy 
can be simply eliminated reducing the angular momentum by the same factor. 

I want to remark that we have used some results of the random Gaussian fields in order 
to calculate the torque (e.g., equation (30)). This could seem a rather strong assumption, 
being us concerned with small scales, where the density field is no more Gaussian. This 
assumption is justified by what previously told, namely by the fact that our calculation of 
the final angular momentum is obtained as an extrapolation of the linear theory and as 
previously quoted this approach gives values of the angular momentum not too different 
from those obtained in numerical simulations. 

As previously quoted, we asssume that from Im on, the ellipsoid has this constant 
angular momentum. Following the procedures 1) and/or 2), we shall be able to get the time 
evolution of the density and that of the collapse velocity. 



In order to apply the model introduced in the previous section to the evolution of the 
collapsing perturbation and solve equations (11)-(14), we need the initial conditions and 
moreover it is necessary to connect this conditions and the time dependence of the shear to 
properties of the initial density field. To begin with, initially, the high density region that 
shall collapse has 6 << 1, and it is contained in a spherical region. Following Eisenstein & 
Loeb (1995), we impose that average mass, density and quadrupole moments of the ellipsoid 
match that of the inner spherical region at the initial time. So defining, as usual, the 
overdensity of the inner region as: 



4. 



Pcirameters, constraints and initial conditions 




(40) 
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the mass of the region is given by: 

Att — Att 

M = —p^R' [I + 8{R)\ = —p^R" [1 + H (41) 

The elhpsoid is chosen to match the previous quantities: it has overdensity 5{R)^ mass, 
M — ^ph [l + 8{R)\ ci;iQ;2«3 and by comparison with equation (41) we get: 

R^ = aia2a2, (42) 

The quadrupole moments, necessary to set g2m,o in equation (35) are obtained from equa- 
tion (27). The initial axes of the elhpsoid are fixed as follows: given a value of 5 (or v) and 
the initial mass, M, we can calculate the radius it! from equation (41). Equation (42), (33), 
(34) make possible to get cti, 0:2, 0:3. The initial density, for the case v — 2, \s 5 — 2 x 10~^, 
and M ~ 2 X IQ^^Mq (since we are concerned with galactic mass scales), and the velocity is 
chosen to be a uniform expansion with the Hubble flow (a pure growing mode). 



The equations of the model described in section 2 were integrated using the Buhrsch- 
Stoer algorithm. We assumed an = 1 universe, a Hubble constant of Hq = 50km/s/Mpc. 
The CDM power spectrum that I adopt is P{k) — AkT'^{k) with the transfer function T{k) 
given in BBKS (equation (G3)): 

T{k) = ^^^Wf^^"^^^ ■ [1 + 3.89g + {W.lqf + {5A6q)' + (6.71)^'/' (43) 

where A is the normalizing constant and q = Here 6 = Per/(l-68p^) represents 

the ratio of the energy density in relativistic particles to that in photons {9 — 1 corresponds 
to photons and three flavors of relativistic neutrinos). The spectrum was smoothed on a 
galactic scale {R ~ O.bh"^ Mpc) and normahzed such that (7(8/i~^Mpc) = 1 at present 
epoch (cTg is the rms value of ^ in a sphere of 8/i~^Mpc). The mass variance present in 
equation (39) can be obtained from the spectrum, P{k), as: 

1 f°° 

a\R) = / dkeP{k)W\kR) (44) 
27r^ Jo 

where W{kR) is a top-hat smoothing function: 

W{kR) = (sin kR - kR cos kR) (45) 

The remaining spectral parameters of equation (39), 7, R^, for the chosen spectrum with the 
above fixed smoothing length, are 7 ~ 0.6, R^ — 0.52. For what concerns the duration of 
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the integration, we followed the suggestions of Eisenstein & Loeb (1995): since the average 
overdensity of the innermost external shell is of the same order of magnitude of that of the 
ellipsoid, the two objects collapse at similar times, or in some cases the inner external shell 
collapses before the long axis of the ellipsoid. To avoid the problem, the integration must be 
stopped at some time before the collapse of the inner tidal shell. This can be accomplished 
(Eisenstein & Loeb 1995) by constraining the initial conditions so that none of the exterior 
shells has an overdensity greater than 95% of the initial density of the eUipsoid. The last 
assumption ensures that the external tidal shells does not collapse before the integration 
ends. As a consequence the inner region behaves as a density peak. We also imposed the 
condition (5(i?sphere) > ^(^1 with i/ > 2, implying that the inner spherical region have high 
overdensity, and finally we follow Bond & Myers (1993a,b), imposing the condition that no 
axis may collapse below 40% of its maximum length, in order to avoid that the dynamics 
approaches the singularity at zero length and to simulate virialization of the corresponding 
axis. 



5. Results 



The results of the calculation involving the evolution of 5 are shown in figures 2-5. 



In figure 2, we plot 5 — in the case of a density peak having height v — S/a{R) — 2. 
The solid line represents the solution of equation (14) in the case L = S = 0. The dashed and 
dotted fines representing the case L = 0, E 7^ and L 7^ 0, S 7^ 0, respectively, were obtained 
using both the procedure 1) and 2) described in section 2. As shown, the shear (see dashed 
line) produces an enhancement of the growth rate of the density of a mass clement. This 
is the growth rate enhancement of the density contrast effect induced by the shear firstly 
recognized by Hoffman (1986a). The shear term, S^, appearing in the equation of evolution 
of the density (equation(14)) is positive definite, so as long as the fiuid is irrotational, the 
growth rate of the density contrast must be enhanced by it, so the effects of the shear are 
present in the linear and nonlinear regime. During the linear phase the ellipsoid expands with 
the universe: the tidal field outside the ellipsoid cancels the effects of the gravitational field 
interior to it. The shear contributes to increase the growth rate of the perturbation. When 
the effect of the angular momentum can no longer be neglected, we see that the situation 
previously described changes (see dotted line). Initially the shear term dominates, but in 
a short time the angular momentum begins to influence the growth of the perturbation, 
counterbalancing the effect of the shear term, E^, and producing a slowing down of the 
growth. As a final result the growth of the density perturbation becomes slower than in the 
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case L — T, — 0. The density contrast at virialization, 5y ~ 60, is reduced with respect 
to the expected value 5^ ~ 178. The value obtained is intermediate between that obtained 
by Peebles (1990) for the half-mass radius, 5^ ~ 30, and that obtained by the modified 
spherical collapse model of Engineer et al. (1998), 6^, ~ 80. We want to remark that in this 
last paper the authors showed that in order to take account of the effects coming from the 
asphericity of a system, one has to add to the equations for the density evolution typical of 
the spherical collapse model an additive term, (1 + 5)(S^ — 2Q^), depending on shear and 
angular momentum of the system, similarly to Nariai & Fujimoto (1972) model. 

The previous result can be interpreted as follows: when the ovcrdcnsity of the ellipsoid 
become considerable , and 6 attains amplitudes of order unity the ellipsoid will begin to 
recoUapse in at least one direction. As shown by Peebles (1980), if we consider the collapse 
of the sphere of equivalent mass, when this reaches the turn-around epoch one of the axis of 
the ellipsoid turns shorter and collapse forming a pancake in which the barions shock and the 
dark matter goes through violent relaxation. In the process the ellipsoid develops nonradial 
motion. The angular momentum L present in equation (11)-(14) becomes not negligible 
and produces the slowing down of density growth shown in the figure. 

In figure 3, we show the same calculation of figure 2, but now we have increased the 
value of the peak height, u = 2.5. As in the previous case, the shear term produces an 
enhancement of the growth rate of the density (dashed line), but this time the effect is 
smaller with respect to that shown in figure 2. This is due to the fact that the shear 
magnitude decreases with increasing peak height (see also Hoffman 1986a, Table 1). As in 
figure 2, the angular momentum acts in the opposite direction to that of shear, but this time 
its effect is reduced (dotted line) with respect to the case u — 2 because of the well known 
L — u anticorrelation effect (Hoffman 1986b). The trend is confirmed by figure 4 representing 
the same calculations of the figures 2-3 but now in the case i/ — 3. 

Summarizing, both L and E reduce their rate of growth with increasing v. rare density 
peaks are in general characterized by a low shear, and then the evolution of the perturbation 
tends to follow the results of the spherical model, when u increases, and as expected also the 
collapse shall be nearly spherical (Bernardeau 1994). 

In order to study the effect of angular momentum and shear on the collapse velocity, we 
calculate the collapse velocity at the epoch of pancaking numerically solving the equations 
of motions for the principal axes of the eUipsoid. Using Barrow & Silk (1981) notation, we 
indicate with XoX{t) and ZoZ{t), the principal axes {xq and Zo are the initial values of the 
axes). We solve equations (11)-(13) to calculate the collapse velocity down the shortest axis 
at the epoch of pancaking in the case of an oblate spheroid (ai = a2 > a^)- This calculation 
is similar to that of Barrow & Silk (1981) with the difference that our approach is a numerical 
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one. Then the collapse velocity at pancaking is: 

^zp = -zjpit) (46) 

(from here on the subscript p indicates that the quantity is evaluated at the pancaking time) . 
The initial conditions are set similarly to section 4 and the equation are solved for several 
values of x^, while = 1. In figure 5, we plot ^ ^""^^ , where Tp = r^Xp is the pancake radius 
(see Barrow & Silk 1981), as function of the ratio of the initial value of the axes, Xo/zo- The 
solid line represents the collapse velocity for an oblate spheroid (ai = 0^2 > 0^3) in the case 
L 7^ 0, E 7^ and u — 3. The dotted and dashed fines plot the result of the same calculation 
but in the case u — 2.5 and u — 2, respectively. The figure shows two trends: 

a) the collapse velocity is reduced with increasing initial flattening (increasing value of Xo)- 
For example for Zo/xo — 0.44 the collapse velocity is reduced to the Hubble velocity in the 
plane of the pancake (HpTp), while in the case of more extreme flattening Zo/xo — 0.125, the 
collapse velocity is reduced by a factor of ^ 6 with respect the previous value. 

b) the collapse velocity is reduced with increasing angular momcntTim acquired by the pro- 
tostructurc. As shown, the collapse velocity is progressively reduced when we go from v = 3 
peaks to u = 2. 

In other words, the slowing down of the rate of growth of density contrast produces a low- 
ering of the peculiar velocity in qualitative and quantitative agreement with Barrow & Silk 
(1981) and Szalay & Silk (1983). 

The result obtained helps to clarify the controversy relative to the previrialization conjec- 
ture. According to this paper, it is surely true that taking account only of the shear, S, 
produces a shortening of the collapse time of non-spherical perturbations, in agreement with 
Hoffman (1986a) and Bertschinger & Jain's collapse theorem. The question is that in the real 
collapse other effects have an important role, namely external tides and the effects of small 
scale substructure. Both Hoffman (1986a) and Bertschinger & Jain (1994) results are vafid 
for a fluid element, which has no substructure by deflnition, while a small scale substructure 
produces a slowing down of the collapse at least in two ways: 

1) encounters between inf ailing clumps and substructure internal to the perturbation (Antonuccio- 
Delogu & Colafrancesco 1994; Del Popolo & Gambera 1997; Del Popolo & Gambera 1999); 

2) tidal interaction of the main proto-structure with substructure external to the perturba- 
tion (Peebles 1990; Del Popolo & Gambera 1998). 

Moreover, it should be pointed out that, as more small-scale power is present the collapse of 
a perturbation may be slowed down in a way that could inhibit the effect of shear. 
Differently from Bertschinger & Jain (1994), our model takes account of the angular momen- 
tum of the system, and then at least of the effects produced by the point 2) quoted above. 
Similarly to Bertschinger & Jain (1994), our model does not take account of the substruc- 
ture internal to the system. This last is a natural limitation of the homogeneous ellipsoid 
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model: a homogeneous ellipsoid cannot represent the substructure of the object. We, how- 
ever, recall that the same shortcoming was present in Peebles (1990) paper: in that paper 
the substructure was suppressed, since it adopted an homogeneous Poisson distribution of 
particles within the protocluster (Peebles 1990). This limit has the effect of underestimating 
the effect of previrialization, and in particular the value of the ovcrdensity at virialization, 
6y (Peebles 1990). In other words, the effects of the slowing down of the collapse obtained 
in this paper (similarly to that of Peebles (1990)) are surely smaller than that we shall find 
if we had used a system having internal substructure, as in the above point 1). 

Before concluding, we want to spend a few words on the impact of the result of the 
paper on our view of structure formation. 

The reduction of the rate of growth of overdensity and collapse velocity has several conse- 
quences on structures formation. To begin with, a first consequence is a change of the mass 
function, the two-point correlation function, and the mass that accretes on density peaks. 
In fact, as several times remarked, the angular momentum acquired by a shell centered 
on a peak in the CDM density distribution is anti-corrclatcd with density: high-density 
peaks acquire less angular momentum than low-density peaks (Hoffman 1986b; R88). A 
greater amount of angular momentum acquired by low-density peaks (with respect to the 
high-density ones) implies that these peaks can more easily resist gravitational collapse and 
consequently it is more difficult for them to form structure. This results in a tendency for 
less dense regions to accrete less mass, with respect to a classical spherical model, inducing 
a biasing of over-dense regions toward higher mass. 

As a result, the number of objects with a < 1 (i.e., large mass) is smaller, since now the col- 
lapse is slowed down, and the mass function is now much below the standard PS prediction 
(Del Popolo & Gambcra 1999; Del Popolo & Gambera 2000; Audit et al. 1997). Even the 
two-point correlation function of galaxies and clusters of galaxies results strongly modified 
(see Del Popolo & Gambera 1999; Del Popolo et al. 1999; Peebles 1993). 

6. Conclusions 

We examined the evolution of non-spherical inhomogeneities in a Einstein-de Sitter uni- 
verse, by numerically solving the equations of motion for the principal axes and the density 
of a dust ellipsoid. We took account of the effect of the mass external to the perturbation 
by calculating the angular momentum transferred to the developing proto-structure by the 
gravitational interaction of the system with the tidal field of the matter becoming concen- 
trated in neighboring proto-structures. 

We showed that for lower values oi u {u — 2) the growth rate enhancement of the density 
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contrast induced by the shear is counterbalanced by the effect of angular momentum acqui- 
sition. For z/ > 3 the effect of angular momentum and shear reduces, and the evolution of 
perturbations tends to follow the behaviour obtained in the spherical collapse model. These 
results corroborate the previrialization conjecture because they show that asphericities and 
tidal torques slow down the collapse of the perturbation after the system detaches from the 
general expansion. 
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Fig. 1. — Comparison of the mean quadrupole moments due to triaxiality (dotted line), 
around a 3cr peak, smoothed on galactic scale, with the sum of the r.m.s. quadrupole 
moments due to the secondary perturbations (solid line). 
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Fig. 2. — Evolution of ellipsoidal density perturbations in an expanding universe as function 
of redshift, for z/ = 2 in the case L — 12 — (solid line), L = 0, E (dashed line) and 
L^O,E^O (dotted line). 




Fig. 3. — Same as figure 2 but now v — 2.5. 




Fig. 4. — Same as figure 2 but now u — 3. 
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Fig. 5. — Collapse velocity for an oblate spheroid {ai — a2 > cks) down the z axis at epoch 
of pancaking (p). Xo and are the initial values of the longest (x) and shortest axis (z), 
Hp and are respectively the Hubble constant and the pancake radius at the pancaking 
epoch. The solid line represents the collapse velocity in the case L 7^ 0, S 7^ 0, and u = 3. 
The dotted and dashed line the same calculation but for i/ = 2.5, = 2, respectively. 



